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ABSTRACT 

We  present  the  results  of  a  new,  comprehensive  investigation  of  the  radio  spectral  index  of  the  Crab 
Nebula  supernova  remnant.  New  data  at  74  MHz  are  combined  with  data  at  327  MHz,  1.5  GHz,  and 
5  GHz.  In  contrast  to  previous  claims,  little  spatial  variation  in  the  spectral  index  is  seen.  In  particular, 
between  327  MHz  and  5  GHz  we  see  no  evidence  of  spectral  steepening  near  the  edge  of  the  Nebula,  the 
“jet,”  or  the  ionized  filaments.  The  rms  limits  on  any  spectral  index  variations  in  these  regions  amount 
to  no  more  than  0.01.  We  believe  that  earlier  reports  of  large  steepening  were  the  result  of  correlator 
bias  and  image  registration  problems.  An  elongated  feature  was  detected  T  northwest  of  the  pulsar, 
which  may  be  a  continuation  of  the  well-known  wisplike  structures  seen  closer  to  the  center  of  the 
Nebula.  At  74  MHz,  we  see  for  the  first  time  evidence  of  free-free  absorption  by  the  thermal  material  in 
the  Crab  Nebula’s  filaments.  Apart  from  some  possible  renewed  acceleration  occurring  in  the  wisps,  the 
dominant  accelerator  of  relativistic  electrons  in  the  Crab  Nebula  is  the  pulsar  itself. 

Subject  headings :  ISM :  individual  (Crab  Nebula)  —  radiation  mechanisms :  nonthermal  — 
radio  continuum :  ISM  —  supernova  remnants 


1.  INTRODUCTION 

The  Crab  Nebula  supernova  remnant  (SNR)  is  one  of  the 
strongest  radio  sources  in  the  sky  and  has  long  been  an 
object  of  intense  scrutiny  by  radio  astronomers.  Neverthe¬ 
less,  a  number  of  controversies  with  important  physical 
implications  continue  to  surround  the  details  of  its  radio 
continuum  spectrum.  While  there  is  general  agreement  that 
the  Crab’s  integrated  spectrum  has  a  slope1 
a.  =  —  0.27  ±  0.04  from  10  MHz  to  the  synchrotron  break 
frequency  near  104  GHz  (Strom  &  Greidanus  1992; 
Kovalenko,  Pynzar’,  &  Udal’tsov  1994),  there  is  consider¬ 
able  dispute  over  the  details  of  its  spatial  distribution 
(Wilson  1972;  Swinbank  1980;  Velusamy,  Roshi,  &  Venu- 
gopal  1992;  Bietenholz  &  Kronberg  1992,  hereafter  BK92). 
Since  the  spectral  index  serves  as  a  useful  measure  of  the 
energy  distribution  of  the  relativistic  electrons  that  generate 
the  synchrotron  emission,  claims  of  spatial  variations  in 
spectral  index  across  the  Nebula  therefore  must  be  inter¬ 
preted  in  terms  of  differences  in  the  underlying  physical 
processes  that  influence  the  energy  of  the  particles. 

The  following  are  the  three  main  areas  where  small  varia¬ 
tions  in  the  radio  spectral  index,  a,  have  been  reported :  In 
the  first  instance,  BK92  report  spectral  index  variations  in 
the  vicinity  of  the  pulsar,  similar  to  the  wisps  found  in  the 
optical  (Scargle  1969;  Hester  et  al.  1995).  Although  they  can 
be  seen  as  weak  arcsecond-width  features  in  the  total  inten¬ 
sity  images,  these  structures  are  most  prominent  in  the  spec- 

1  We  define  the  radio  spectral  index,  a,  so  that  Sv  oc  v“,  where  Sv  is  the 
flux  density  at  frequency  v. 


tral  index  images.  They  have  a  flatter  spectral  index  than 
the  Nebula,  which  suggests  that  they  are  regions  of  renewed 
particle  acceleration.  These  spectral  index  ridges  are  located 
in  the  same  region  but  are  not  positionally  coincident  with 
the  optical  wisps.  They  have  the  same  concave,  arclike 
appearance  as  the  optical  wisps  but  are  somewhat  longer. 
These  ridges  were  interpreted  as  manifestations  of  the  wind 
shock  in  the  Crab  in  BK92,  and  Kronberg  et  al.  (1993) 
modeled  them  as  such  with  a  tilted  ring  or  torus  geometry. 

Second,  Velusamy  et  al.  (1992)  have  claimed  that  the 
spectrum  of  the  filaments  in  the  body  of  the  Nebula  is  differ¬ 
ent  than  that  of  the  intrafilament  plasma  and  that  there  is  a 
trend  of  steepening  filament  spectral  index  as  a  function  of 
radius  in  the  Nebula.  No  such  variations  were  reported  in 
BK92,  however.  If  true,  these  spectral  variations  must  be 
reconciled  with  interpretations  of  the  physical  nature  of  the 
filaments  and  the  source  of  their  radio  emission  (Reynolds 
1988).  This  would  imply  a  different,  or  at  least  an  addi¬ 
tional,  acceleration  process  for  the  relativistic  electrons  in 
the  filaments. 

Third,  several  groups  have  reported  a  steepening  of  the 
spectral  index  near  the  edge  of  the  Nebula  (Velusamy  & 
Sarma  1977;  Agafonov  et  al.  1987).  These  measurements, 
however,  are  difficult  to  make,  and  other  groups  report  no 
steepening  of  os  toward  the  edges  of  the  Nebula  (Wilson 
1972;  Swinbank  1980;  Trushkin  1986;  BK92).  If  such  a 
steepening  is  real,  then  it  could  be  interpreted  as  evidence 
for  a  shell  of  emission  around  the  Crab,  analogous  to  the 
steeper  spectrum  shell-type  emission  typically  observed 
from  classical  shell-type  SNRs  (e.g.,  Cas  A).  The  question  of 
whether  or  not  such  a  shell  exists  superposed  on  the  plerion 
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TABLE  1 

Details  of  the  Observing  Sessions 


Frequencies  VLA 

(MHz)  Epoch  Arrays  Used  Reference 


4885,4625 .  1987  BCD  Bietenholz  &  Kronberg  1991 

1515,1410 .  1987  ABCD  Bietenholz  &  Kronberg  1991 

327 .  1992  BC  Frail  et  al.  1995 

74  .  1993  ABC  Kassim  et  al.  1993 


is  extremely  important,  since  if  it  does,  this  would  suc¬ 
cessfully  bring  to  an  end  the  search  for  the  classical  radio 
shell  exterior  to  the  presently  observed  Crab  SNR.  The  lack 
of  such  a  shell  has  led  previous  workers  to  important  con¬ 
clusions  regarding  the  energy  of  the  supernova  explosion, 
the  efficiency  of  its  blast  wave  in  accelerating  relativistic 
particles,  and  the  nature  of  the  circumstellar  environment 
surrounding  the  SN  precursor  star  (Nomoto  1985,  1987; 
Frail  et  al.  1995,  hereafter  F95).  If  the  existence  of  a  shell 
were  confirmed,  all  conclusions  drawn  from  those  negative 
searches  would  be  in  error. 

In  order  to  tackle  these  controversies,  we  bring  together 
in  this  paper  a  comprehensive  body  of  observational 
material  and  present  important  new  subarcminute- 
resolution  74  MHz  observations.  This  radio  database, 
encompassing  an  almost  two  decade  range  in  frequency, 
allows  us  to  conduct  the  most  sensitive  study  of  the  magni¬ 
tude  and  spatial  distribution  of  the  Crab’s  radio  spectral 
index  conducted  to  date.  We  then  use  this  unique  body  of 
evidence  to  address  the  observational  controversies  out¬ 
lined  above  and  to  address  the  physical  implications  of 
these  findings. 

2.  OBSERVATIONS  AND  DATA  REDUCTION 

We  used  several  sets  of  radio  observations  as  described  in 
Table  1.  All  data  were  taken  at  the  NRAO  Very  Large 
Array  (VLA).2  These  include  74  MHz  observations  that 
currently  utilize  only  eight  VLA  antennas.  The  74  MHz 
data  comprised  snapshot  observations  from  three  A,  two  B, 
and  one  C  array  configuration  spanning  several  years  in 
order  to  obtain  sufficient  u-v  coverage.  The  nonstandard 
procedures  employed  to  reduce  these  data  are  described  in 
Kassim  et  al.  (1993).  The  image  obtained  at  74  MHz  is 
shown  in  Figure  1. 

Several  problems,  which  have  not  been  addressed  fully  in 
the  past,  complicate  the  determination  of  the  spectral  index 
from  interferometer  data  and  especially  from  the  VLA. 

First,  the  images  at  different  frequencies  need  to  be  regis¬ 
tered  accurately.  Most  radio  images  have  had  the  antenna 
phases  improved  by  self-calibration  (see,  e.g.,  Pearson  & 
Readhead  1984;  Perley  1989).  Once  this  has  been  done,  the 
absolute  position  information  is  degraded.  To  make  a  spec¬ 
tral  index  map,  however,  the  two  maps  must  be  registered 
accurately,  and  so  a  way  of  recovering  the  absolute  position 
is  needed.  Typically  this  problem  can  be  ignored,  for  if  the 
starting  model  is  reasonable,  the  shift  in  position  is  small, 
and  registering  maps  by  the  nominal  positions  is  sufficient. 
Nonetheless,  shifts  of  up  to  one-fourth  of  a  beamwidth  can 
occur.  In  the  case  of  a  completely  resolved  source,  a  shift  of 


2  The  NRAO  VLA  is  a  facility  of  the  National  Science  Foundation 
operated  under  cooperative  agreement  by  Associated  Universities,  Inc. 


only  1/10  of  a  beam  can  change  the  map  amplitude  at  a 
particular  pixel  position  by  up  to  25%,  which  can  cause 
substantial  spurious  spectral  index  variations,  especially 
near  bright  gradients  such  as  those  near  the  edges  of  SNRs. 

Images  of  the  Crab  are  complicated,  and  registering  them 
to  less  than  1/10  beamwidth  is  difficult  in  the  presence  of 
noise.  One  means  of  accomplishing  this  is  by  using  the 
MIRIAD  program  IMDIFF  (Dickel  et  al.  1988),  which  uses 
the  maximum  entropy  method  of  Tan  &  Gull  (1985)  to  find 
the  translation  required  to  make  the  first  image  resemble 
the  second  one  as  closely  as  possible  in  a  least-squares  sense 
(see  Bietenholz  et  al.  1991  for  an  application  of  this  method 
to  radio  images  of  the  Crab).  The  other  approach  is  to 
self-calibrate  the  data  in  phase  using  a  model  of  known 
positional  offset,  which  will  serve  to  bring  the  positions  to 
the  reference  frame  defined  by  the  model  (i.e.,  using  the 
image  at  one  frequency  to  phase  calibrate  the  data  at 
another  frequency). 

A  much  less  tractable  problem,  peculiar  to  the  VLA,  is 
the  correlator  nonlinearity  (Kulkarni  &  Heiles  1980).  The 
correction  needed  to  compensate  for  this  is  known  as  the 
quantization  or  van  Vleck  correction,  and  it  is  not  done  at 
the  VLA.  To  our  knowledge,  none  of  the  previous  VLA 
results  have  addressed  this  problem.  The  output  of  the 
digital  correlator  employed  at  the  VLA  is  a  known  function 
of  the  true  analog  correlation  factor.  For  small  correlated 
flux  densities,  this  function  is  very  nearly  linear,  and  this 
linear  approximation  is  what  is  used  at  the  VLA  (D’Addario 
1989).  However,  for  objects  as  strong  as  the  Crab,  the  rela¬ 
tion  between  the  digital  correlation  and  the  true  analog 
correlation  is  no  longer  linear,  and  the  derived  correlated 
flux  densities  will  be  in  error.  This  error  increases  with  the 
observed  flux  density.  It  affects  the  real  and  imaginary  part 
of  the  correlation  coefficient  independently,  and  since  the 
uncalibrated  phase,  i.e.,  that  of  the  correlated  flux  density,  is 
random,  it  will  cause  random  phase  and  amplitude  errors. 
The  observed  amplitudes,  however,  will  always  be  lower 
than  the  true  ones,  so  the  resulting  maps  will  be  biased.  For 
our  data,  at  5  GHz,  these  errors  can  be  as  high  as  4%  in 
amplitude. 

Unfortunately,  the  full  information  needed  to  correct  for 
this  effect  is  not  as  yet  preserved  at  the  VLA,  so  it  is  difficult 
to  accurately  recover  the  true  correlation  after  the  fact.  This 
is  less  likely  to  affect  our  images  at  74  MHz  because  the 
correlator  nonlinearity  is  less  severe  at  this  frequency,  and 
our  signal-to-noise  ratio  is  considerably  lower.  For  the  data 
at  other  bands,  we  investigated  the  effects  of  this  problem  by 
making  two  sets  of  spectral  index  maps.  The  first  is  the 
“standard”  version  that  is  made  in  the  straightforward 
way,  with  the  maps  being  registered  by  IMDIFF.  Second, 
we  made  “  careful  ”  spectral  index  maps,  in  which  we  decon¬ 
volved  using  only  data  at  u-v  distances  >  1.6  k/.,  where  the 
correlated  flux  density  is  low  enough  that  the  error  in 
neglecting  the  correlator  bias  is  minor.  The  registration  was 
achieved  by  also  self-calibrating  in  phase  using  a  common 
model.  We  describe  the  production  of  these  images  in  more 
detail  below. 

Whether  any  significant  bias  is  introduced  into  the  maps 
by  the  correlator  nonlinearity  is  difficult  to  estimate.  The 
process  of  deconvolution  routinely  recovers  the  flux  density 
in  the  center  of  the  u-v  plane  starting  from  an  estimate  of 
zero,  so  one  might  think  that  it  should  have  no  problem 
recovering  a  few  percent  of  missing  flux  density.  On  the 
other  hand  the  deconvolution  process  is  completely  free  to 
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Fig.  1. — The  74  MHz  image  of  the  Crab,  with  the  pulsar  subtracted.  The  FWHM  of  the  restoring  beam  was  25".  The  peak  flux  density  is  26  Jy  beam-1, 
and  the  noise  level  is  0.5  Jy  beam- '  .  The  contours  are  drawn  at  11.3%,  16%,  22.6%,  32%,  and  45.3%  of  the  peak  flux  density. 


adjust  the  flux  density  in  the  center  of  the  u-v  plane,  whereas 
in  regions  where  there  are  measurements,  deconvolution 
will  attempt  to  fit  the  biased  measurements.  We  found  that 
the  apparent  noise  levels  in  the  maps  deconvolved  by  our 
“  careful  ”  process  was  lower  than  the  “  standard  ”  ones.  This 
suggests  that  there  is  discernible  noise  introduced  by  the 
missing  correction,  since  one  does  not  expect  the  rms  to 
decrease  when  one  reduces  the  amount  of  data.  This  would 
suggest  that  the  bias  might  also  be  significant. 

We  will  now  describe  more  fully  the  production  of  the 
careful  spectral  index  images:  Because  the  measurements 
near  the  center  of  the  u-v  plane  contain  the  most  correlated 
flux,  they  are  most  severely  affected  by  the  correlator  non¬ 
linearity.  We  made  maps  without  using  the  data  that  were 
compromised  in  this  way.  We  chose  to  use  a  u-v  distance 
cutoff  rather  than  a  flux  density  cutoff  because  its  effect  on 
the  images  is  easier  to  understand.  The  resulting  spectral 
index  images  will  contain  little  information  on  angular 
scales  larger  than  217,  i.e.,  half  the  size  of  the  Nebula.  At  5 
GHz,  which  is  the  worst  case,  the  correlated  flux  density  at 


this  u-v  distance  is  ~20  Jy,  and  the  remaining  correlator 
nonlinearity  errors  will  be  <  0.05%. 

Since  the  Crab  Nebula  is  considerably  larger  than  the 
largest  scale  of  217  that  we  are  sampling,  we  will  recover  the 
large-scale  structure  by  deconvolving  with  a  default 
(support)  in  our  maximum  entropy  deconvolutions 
(Astronomical  Image  Processing  System  [AIPS]  task 
YTESS;  see  Cornwell  &  Evans  1985;  Cornwell  1988).  We 
will  use  the  same  default  (appropriately  scaled)  for  all  of  our 
observing  frequencies,  which  is  justified  because  the  spectral 
index  is  quite  uniform  over  the  Nebula  (Wilson  1972; 
BK92).  We  will  concentrate  on  mapping  the  difference  to 
the  Nebula’s  average  spectral  index  rather  than  its  absolute 
spectral  index,  which  can  be  more  reliably  determined  using 
other  methods:  Kovalenko  et  al.  (1994)  have  recently  deter¬ 
mined  it  to  be  a  =  —  0.27  ±  0.04.  Using  the  same  default 
map  will  bias  our  maps  toward  looking  similar,  and  thus 
any  spatial  variation  in  spectral  index  will  be  required  by  the 
data  (this  is  the  “  forced  ”  YTESS  technique  of  Anderson  & 
Rudnick  1993). 
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Fig.  2. — The  relative  radio  spectral  index  between  327  MHz  and  1.5  GHz  (Aa^J)  in  units  of  0.001.  These  maps  show  the  difference  of  the  spectral  index  to 
the  mean  spectral  index  of  the  Nebula.  The  FWHM  of  restoring  beam  is  16".  The  contours  are  the  0.5%,  1%,  and  2%  contours  in  total  intensity  of  the  327 
MHz  image,  (a)  The  careful  version,  made  by  deconvolving  with  a  default  and  using  only  data  at  u-v  distances  >  1.6  kX.  It  contains  little  information  on 
angular  scales  >2'J  but  is  not  corrupted  by  correlator  nonlinearity.  ( b )  The  standard  version,  which  was  made  without  a  default  and  contains  information  on 
all  angular  scales  but  may  be  slightly  biased. 


The  default  map  was  produced  as  follows:  At  each  fre¬ 
quency  (except  at  74  MHz)  we  mapped  the  data  without 
any  u-v  restrictions.  The  resulting  restored  maps  (i.e.,  con¬ 
volved  with  an  appropriate  “clean”  beam  and  with  the 
residuals  added)  were  corrected  for  the  primary  beam 
pattern  and  for  the  expansion  of  the  Crab  (0.13%  yr_1; 
Bietenholz  et  al.  1991).  They  were  then  convolved  to  a 
common  resolution  of  16"  x  16"  and  flux-scaled  to  a 
common  frequency  using  the  Baars  et  al.  (1977)  spectral 
index  of  a  =  —0.3.  The  recent  Kovalenko  et  al.  (1994)  value 
of  a  =  —  0.27  is  consistent  with  the  older  Baars  et  al.  value 
to  within  the  uncertainty,  so  using  the  latter  should  not 
affect  our  results.  These  maps  were  then  averaged,  and  the 
result  was  used  as  the  default  map  in  all  the  subsequent 
deconvolution  processes.  For  each  deconvolution,  it  was 
again  scaled  back  to  the  flux  density  and  epoch  appropriate 
for  that  frequency. 


Furthermore,  in  order  to  ensure  accurately  registered 
maps,  this  same  (scaled  and  expanded)  default  image  was 
used  to  self-calibrate,  in  phase  only,  the  u-v  data  prior  to 
deconvolution.  This  process  has  the  side  effect  of  further 
biasing  the  maps  toward  being  as  similar  as  is  allowed  by 
the  data.  In  other  words,  it  will  minimize  any  spatial  varia¬ 
tions  in  the  resulting  spectral  index  map.  However,  our 
dynamic  range  is  high,  and  we  will  show  that  this  bias  is 
negligible  and  that  our  careful  spectral  index  maps  can  be 
used  with  confidence. 

Once  self-calibrated  in  this  way,  the  u-v  data  were  trans¬ 
formed  into  the  image  plane,  using  only  data  at  u-v  dis¬ 
tances  greater  than  1.6  kl  and  then  deconvolved  using  the 
AIPS  task  VTESS  and  the  appropriately  scaled  default 
map. 

The  resulting  restored  maps  were  convolved  to  a 
common  beam  size,  contracted  to  a  epoch  1987.4,  and  then 
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spectral  index  maps  were  made  in  the  usual  way.  The  uncer¬ 
tainties  were  derived  from  empty  regions  of  the  maps  near 
the  Nebula  (especially  at  5  GHz;  the  noise  in  a  map  cor¬ 
rected  for  primary  beam  pattern  increases  with  distance 
from  the  phase  center).  These  noise  estimates  were  consis¬ 
tent  with  the  differences  between  4625  and  4885  MHz  maps. 
At  74  MHz  the  pulsar  is  quite  strong,  and  we  first  sub¬ 
tracted  it  from  the  fully  calibrated  u-v  data.  Furthermore, 
the  74  MHz  observations  have  much  lower  signal-to-noise 
ratios  and  fewer  antennas,  and  the  correlator  nonlinearity  is 
less  severe,  so  we  did  not  make  a  careful  spectral  index  map 
involving  this  frequency. 

3.  RESULTS 

The  relative  spectral  index  between  327  MHz  and  1.5 
GHz  (AofoJ)  is  shown  in  Figure  2,  while  Figure  3  shows  the 
relative  spectral  index  between  327  MHz  and  4.9  GHz 
(Aojo  3),  both  at  a  resolution  of  16"  x  16".  These  are  spectral 
indices  relative  to  the  mean  spectral  index  of  the  Nebula 


(—0.27);  we  are  not  sensitive  to  the  total  flux  densities,  so 
our  integrated  spectral  indexes  are  not  constrained  by  our 
data.  Our  reduction  process  will  make  them  very  close  to 
the  value  of  —  0.3,  which  we  assumed  for  the  default  maps. 
In  addition,  Figure  4  shows  two  profiles  through  the  AojqJ 
spectral  index  map.  The  profiles  go  through  the  pulsar  posi¬ 
tion  at  position  angles  —  45°  and  +  45°. 

We  have  included  both  the  careful  and  the  standard  maps 
in  order  to  assess  the  significance  of  the  bias  (see  §  2)  and  the 
degree  of  success  we  have  had  in  treating  it.  Both  maps 
reveal  that  the  spectral  index  of  the  Nebula  is  quite  uniform. 
The  rms  over  the  whole  Nebula  of  AoeqJ,  for  example,  is 
only  0.01  in  a  when  one  excludes  the  region  near  the  edge  of 
the  Nebula  where  the  signal-to-noise  ratio  is  low.  The  differ¬ 
ences  between  the  careful  and  the  standard  maps  are  negli¬ 
gible  in  the  case  of  AojqJ ;  the  rms  of  the  difference  in  a  is 
0.008.  There  is  no  evidence  of  spectral  steepening  of  any 
significant  magnitude  near  the  edges  of  the  Nebula.  An 
“  apparent  ”  steepening  is  visible  at  some  locations,  near  the 
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Fig.  3  a 

Fig.  3. — Same  as  Fig.  2,  but  for  the  relative  radio  spectral  index  between  327  MHz  and  5  GHz(Aa„  3)  in  units  of  0.001  (see  text,§  2) 


nebular  boundaries.  However,  these  occur  at  brightness 
levels  less  than  1%  of  the  peak  brightness  of  the  Crab 
Nebula  (at  327  MHz),  where  the  uncertainty  in  Aa  is  large. 

The  pulsar  is  the  strongest  feature  in  these  spectral  index 
maps.  It  appears  as  a  point  source  of  steeper  spectral  index 
visible  in  the  AocqJ  and  Aajj.3  maps.  It  shows  up  as  a 
steepening  of  only  0.07  on  the  Aajj.3  map.  This  translates 
into  a  pulsar  flux  density  of  1.9  Jy  at  327  MHz,  which  is 
consistent  with  the  value  of  1.3  Jy  given  by  Lorimer  et  al. 
(1995),  since  pulsar  fluxes  are  known  to  be  quite  variable. 
Since  the  spectral  index  at  the  pulsar  did  not  change  notice¬ 
ably  between  the  careful  and  standard  maps,  any  bias  intro¬ 
duced  by  our  careful  procedure  is  small.  All  the  maps  also 
show  a  slight  flattening  in  a  region  of  ~1'  to  2'  in  size 
surrounding  the  pulsar  and  extending  slightly  farther  to  the 
northwest.  The  magnitude  of  this  effect  ranges  from 
Aa  ~  0.05  on  the  standard  Aajj.3  maP  to  Aa  —  0-02  on  the 
careful  Aajj.3  and  the  AajjJ  maps.  It  is  notably  larger  on  the 
standardAa;)  3  map  than  on  the  careful  one.  Since  the 
angular  scales  involved  are  those  that  could  be  biased  by 
the  correlator  nonlinearity,  which  is  worst  at  5  GHz,  we 


suspect  that  the  estimate  of  Aa  =  0.02  is  more  reliable.  This 
effect  is  gradual,  and  remaining  systematic  sources  of  error 
are  probably  at  about  this  level.  We  note  that  the  correlator 
nonlinearity  would  tend  to  reduce  the  observed  flux  density 
especially  at  higher  frequencies  and  at  large  angular  scales, 
which  could  cause  a  slight  spurious  spectral  index 
steepening  to  appear  toward  the  edges  of  the  Nebula  in 
VLA  images. 

We  emphasize  that  our  conclusions  given  above  are 
largely  independent  of  whether  we  examine  the  careful  or 
the  standard  spectral  index  maps;  this  suggests  that  the  bias 
introduced  by  the  correlator  nonlinearity  is  in  fact  small. 
Our  relatively  uniform  spectral  index  images  suggest  that 
the  visibilities  at  two  frequencies  should  differ  only  by  a 
scale  factor.  We  examined  the  gridded  visibility  data,  and 
they  show  the  expected  linear  relationship  between  327  and 
1.5  GHz  at  u-v  distances  >  1.6  k/..  The  formal  scatter  is  ~2 
cr,  which  is  compatible  with  the  small  spectral  index  differ¬ 
ences  found  (since  the  above  process  accounted  for  neither 
residual  calibration  errors  nor  the  primary  beam  pattern, 
we  expect  the  scatter  about  a  linear  relationship  to  be  some- 
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Fig.  3  b 


what  larger  than  1  even  in  the  case  of  a  completely  uniform 
spectral  index). 

The  “wisps”  that  were  visible  in  the  higher  resolution 
spectral  index  map  of  BK82  are  barely  discernible  at  this 
resolution.  This  is  simply  because  our  beam  area  is  ~70 
times  that  of  BK82.  Another  feature  is  consistently  visible 
on  all  of  the  Aosq  3  and  Aa£;§  maps:  it  is  a  lane  of  steeper 
spectral  index  ~1'  to  the  northwest  of  the  pulsar.  This 
feature  has  an  AaJ;§  that  is  ~0.03  steeper  than  the  sur¬ 
rounding  region. 

The  spectral  index  of  the  Crab  Nebula’s  “jet  ”  is  insignifi¬ 
cantly  steeper  than  that  of  the  body  of  the  Nebula.  For¬ 
mally,  the  AaJ;§  of  the  jet  is  +  0.05  ±  0.07  relative  to  the  rest 
of  the  Nebula.  We  conclude  that  the  jet  seems  to  have  the 
same  spectral  index  as  the  rest  of  the  Nebula.  This  is  consis¬ 
tent  with  the  star  trail  model  of  Cox,  Gull,  &  Greene  (1991), 
which  explains  the  jet  as  the  star  trail  of  the  Crab’s  progeni¬ 
tor.  The  star  trail  has  cooled  into  a  hollow  tube,  and  the 
lower  pressure  inside  the  tube  causes  the  relativistic  plasma 
from  the  Nebula  to  expand  more  rapidly  there. 


Figure  5  shows  the  map  of  the  spectral  index  between  74 
MHz  and  1.5  GHz  ((x^m),  this  time  at  the  slightly  coarser 
resolution  of  27"  x  25".  Because  of  the  lower  signal-to-noise 
ratio  of  the  74  MHz  map  and  because  the  correlator  nonlin¬ 
earity  is  small  at  this  lowest  frequency,  we  expect  the  errors 
caused  by  not  correcting  for  it  will  be  less  than  the  noise. 
The  pulsar  is,  of  course,  not  visible  in  the  74  MHz  map  or  in 
the  0C74M  map  because  it  has  already  been  removed  from  the 
74  MHz  data.  However,  some  significant  structure  is  visible 
in  this  map.  The  most  distinct  feature  on  this  map  is  at 
R.A.  =  05h31m27!7,  decl.  =  21°58'42':0  and  has  a  spectral 
index  that  is  ~  0.3  flatter  than  the  rest  of  the  Nebula.  It  is 
most  likely  a  low-frequency  absorption  feature  because  no 
corresponding  feature  is  visible  at  frequencies  above  327 
MHz.  Its  reality  was  tested  by  deconvolving  the  74  MHz 
data  using  an  appropriately  scaled,  translated,  and  con¬ 
volved  version  of  the  1.5  GHz  map  as  a  default — once  again 
biasing  the  maps  to  be  as  similar  as  is  allowed  by  the  data — 
and  it  did  not  diminish  in  strength.  Furthermore,  since  the 
74  MHz  data  is  derived  from  a  number  of  observing  runs 
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Fig.  4. — Two  profiles  of  the  A aJJ  spectral  index  map  (careful  version). 
The  profiles  go  through  the  pulsar  position  at  the  indicated  position  angles 
(  —  45°  and  +45°).  Squares  indicate  AaJ;3,  while  the  thin  line  indicates  the 
total  intensity  at  327  MHz,  included  for  comparison. 

and  the  feature  in  question  is  not  confined  to  only  one  of 
them,  it  is  almost  certainly  neither  an  instrumental  nor  a 
deconvolution  effect. 

4.  DISCUSSION 

We  see  no  steepening  of  the  spectral  index  near  the  edge 
of  the  Nebula,  which  confirms  the  result  presented  in  BK92. 
We  find  no  support  for  the  hypothesis  proposed  by 
Yelusamy  &  Sarma  (1977)  that  the  emission  near  the  edge 
of  the  Nebula  comes  from  a  supernova  shell  that  is  super¬ 
posed  on  the  plerionic  emission  from  the  rest  of  the  Nebula. 
On  our  very  sensitive  spectral  index  maps  between  74  MHz 
and  5  GHz,  we  find  that  the  mean  spectral  index  near  the 
edge  of  the  Nebula3  is  —  0.01  +  0.04  flatter  than  that  in  the 
interior  of  the  Nebula.  The  spectral  index  between  74  MHz 
and  1.5  GHz  (a)45M)  shows  no  steepening  to  the  edge  of  the 
Nebula  either,  so  the  spectrum  near  the  edge  maintains  its 
power  law  form  at  least  down  to  74  MHz.  The  uniformity  of 
the  spectral  index  throughout  the  Nebula  suggests  that  all 
the  synchrotron-emitting  electrons  have  the  same  origin, 
which  implies  that  they  are  all  accelerated  by  the  pulsar. 

Velusamy  et  al.  (1992)  claim  that  the  spectral  index  of  the 
filaments  ( af ),  as  opposed  to  that  of  the  smooth  continuum 
emission,  steepens  gradually  toward  the  edge  of  the  Nwith 
uf  increasing  smoothly  by  0.3  from  the  center  to  the  edge  of 

3  The  region  “  near  the  edge  ”  was  defined  as  the  region  between  the  2% 
and  20%  contours  of  the  5  GHz  total  intensity  map,  on  which  the  peak  flux 
density  was  5.2  Jy  beam-1.  The  quoted  uncertainties  are  the  rms  of  a  per 
beamwidth  over  those  regions. 


the  Nebula.  Such  a  steepening  would  be  apparent  in  our 
data  even  though  we  do  not  separate  the  filament  emission 
from  the  continuum;  Yelusamy  et  al.  (1992)  find  that  the 
filaments  account  for  20%  of  the  total  flux  density  at  327 
MHz,  and  since  their  filling  factor  is  <  we  estimate  that 
their  contribution  to  the  surface  brightness  at  filament  loca¬ 
tions  should  be  >  40%.  This  would  cause  the  a  of  the  com¬ 
bined  emission  from  the  continuum  and  filaments  to  change 
by  at  least  0.1,  which  is  considerably  larger  than  any  spec¬ 
tral  index  change  we  observe. 

A  very  slight  and  gradual  steepening  of  the  spectrum  is 
observed  as  one  moves  away  from  the  pulsar.  The  magni¬ 
tude  of  this  effect  is  ~  0.02  in  a.  While  the  reality  of  this 
feature  was  questioned  in  §  3,  we  may  expect  some 
steepening  to  occur,  since  the  radio-emitting  electrons  are 
moving  outward  and  losing  energy  by  synchrotron  radi¬ 
ation.  However,  unlike  in  much  longer  lived  extragalactic 
sources,  where  synchrotron  aging  can  dominate  (Carilli  et 
al.  1991),  the  magnitude  of  this  effect  may  be  too  small  to 
account  for  this  slight  observed  steepening  of  0.02  in  a.  If  we 
assume  an  (average)  nebular  field  of  5  x  10“ 5  G  and 
that  the  electrons  radiate  all  their  power  at  vm  = 
1.22  (B/ 1  G-1)y2  MHz,  we  can  estimate  the  change  in  spec¬ 
tral  index  of  only  <  0.01  over  the  lifetime  of  the  Nebula. 

The  visible  Crab  Nebula  is  an  expanding  bubble  of 
thermal  filaments  and  relativistic  gas.  It  must,  at  some  level, 
interact  with  the  exterior  medium  into  which  it  expands. 
This  interaction  would  probably  take  the  form  of  a  shell.  If 
it  is  expanding  into  the  ISM,  this  shell  would  presumably  be 
similar  to  other,  shell-like  supernova  remnants.  At  327 
MHz,  F95  reported  no  extended  shell  emission  exterior  to 
the  Nebula  and  thus  no  evidence  for  a  fast  hydrogen 
envelope  around  the  Crab.  If  there  is  no  fast  envelope 
around  the  visible  Nebula,  however,  then  the  interface 
between  the  Nebula  and  the  exterior  medium  is  at  the 
observed  edge  of  the  Nebula.  This  interface  presumably 
involves  a  shock  that  could  accelerate  particles.  The  emis¬ 
sion  from  this  interaction  would  thus  be  superposed  on  the 
visible  Nebula.  Shock  acceleration  typically  produces  radio 
spectra  steeper  (a  ~  —  0.5)  than  that  of  the  pulsar- 
accelerated  interior  of  the  Nebula,  so  if  there  were  a  signifi¬ 
cant  quantity  of  shock-accelerated  electrons  near  the  outer 
boundary,  a  steepening  of  the  radio  spectral  index  would  be 
expected.  We  find  no  such  steepening  or  any  “rim”  around 
the  outer  edges  of  the  Nebula. 

Taking  the  limit  of  the  change  in  a  near  the  edge  to  be 
0.03,  we  find  that  the  maximum  residual  surface  brightness 
with  a  =  —  0.5  near  edge  of  the  Nebula  is  <  0.15  Jy  beam- 1 
at  327  MHz.  This  would  represent  <  60  Jy  if  spread  evenly 
over  the  Nebula,  which  would  entail  a  fraction  of  less  than 
2%  of  the  total  radio  luminosity  of  the  Crab  (integrated 
from  107  to  1011  Hz)  and  even  less  if  the  a  =  —0.5  emission 
displayed  a  shell  morphology.  This  would  correspond  to  a 
minimum  relativistic  particle  energy  in  the  putative  shell  of 
less  than  few  times  1048  ergs  (using  standard  assumptions 
from  Pacholczyk  1970).  While  shell-type  supernova  rem¬ 
nants  typically  convert  between  1%  and  10%  of  the  blast 
energy  into  shock-accelerated  synchrotron-emitting  par¬ 
ticles  (Duric  et  al.  1995),  the  Crab’s  blast  wave  is  expected  to 
have  converted  a  far  smaller  fraction  of  its  energy  into  par¬ 
ticles  for  two  reasons :  (1)  The  Crab  is  not  yet  in  the  Sedov 
phase  of  its  evolution,  and  calculations  show  the  efficiency 
of  particle  acceleration  is  very  low  in  the  free-expansion 
phase  (Drury,  Markiewicz,  &  Volk  1989;  Markiewicz, 
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Fig.  5. — The  radio  spectral  index  between  74  MHz  and  1.5  GHz  (a,^,)  in  units  of  0.001.  The  size  of  the  restoring  beam  was  27"  x  25"  elongated  at 
position  angle  —  33°.  The  contours  are  in  total  intensity  at  1.5  GHz,  drawn  at  0.5%,  1%,  2%,  4%,  and  8%  of  the  peak  flux  density. 


Drury,  &  Volk  1990).  In  any  case,  since  the  shell  is  observed 
to  be  accelerating  (Bietenholz  et  al.  1991),  it  cannot  have 
converted  a  great  deal  of  its  original  energy  into  relativistic 
particles.  (2)  The  Crab  supernova  may  well  have  been  an 
anomalously  low-energy  event  (Nomoto  1985,  1987;  Pols 
1994;  F95).  In  the  case  that  the  visible  Crab  is  expanding 
into  a  halo  of  pre-SN  wind  or  ejecta,  as  hypothesized  by 
Hester  et  al.  (1996),  we  would  also  expect  <  1%  of  the 
energy  of  the  expanding  bubble  to  have  been  converted  in 
accelerated  particles,  since  the  relative  velocity  of  the  syn¬ 
chrotron  bubble  is  lower  relative  to  the  exterior  medium. 
Again,  we  would  then  not  expect  to  see  synchrotron  emis¬ 
sion  from  the  shock-accelerated  electrons  at  the  edge  of  the 
synchrotron  bubble. 

The  lane  of  steeper  spectral  index  mentioned  in  §  2  is  ~  V 
to  the  northwest  of  the  pulsar.  This  corresponds  in  position 
to  the  moving  feature  found  in  the  radio  (BK92),  and  to  the 
moving  “dark  lane”  visible  in  the  optical  (Oort  &  Walra- 
ven  1956;  Scargle  1969).  This  feature  is  moving  outward, 
and  in  the  radio,  it  seems  to  have  a  sharp  outward  edge 
(BK92).  Since  our  327  MHz  data  was  taken  several  years 
after  the  data  at  1.5  and  5  GHz,  we  would  expect  such  a 


moving  feature  to  generate  a  spurious  steeper  spectral  index 
feature.  We  can  roughly  estimate  the  feature’s  proper 
motion  to  be  ~0"4  +  0"2  yr-1,  representing  a  speed  of 
~  4000  km  s-1.  This  is  marginally  consistent  with  the 
optical  proper  motion  of  0"7±0''l  yr-1  (Scargle  1969). 
Possibly,  the  feature  is  slowing  down,  since  BK92  find  a 
marginally  slower  speed  (0"6  yr-1)  in  1985.  This  is  consis¬ 
tent  with  the  feature  being  embedded  in  the  large-scale  flow 
in  the  Nebula.  Higher  resolution  observations  would  be 
needed  to  accurately  determine  the  speed  of  the  feature. 

5.  LOW-FREQUENCY  ABSORPTION 

As  discussed  above,  the  significant  feature  southwest  of 
the  pulsar  on  the  spectral  index  map  has  an  index  that 
is  0.3  flatter  than  the  rest  of  the  Nebula  and  that  suggests 
there  is  low-frequency  absorption  toward  this  location.  This 
feature  has  not  been  seen  in  prior  spectral  index  maps,  and 
it  is  the  availability  of  a  subarcminute-resolution  image  at 
74  MHz  that  has  enabled  its  detection.  Comparison  with 
the  three-dimensional  [O  m]  data  of  Lawrence  et  al.  (1995) 
indicates  that  it  corresponds  to  the  location  of  the  brightest 
optical  knot  with  negative  radial  velocity.  This  places  the 
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Fig.  6. — Plot  of  intensity  of  [O  hi]  25007  A  emission  with  negative 
radial  velocity,  i.e.,  that  originating  from  the  side  of  the  Nebula  toward  us 
vs.  the  74  MHz  to  1.5  GHz  spectral  index  (a,i5M).  Optical  data  kindly 
supplied  by  S.  Lawrence  (Lawrence  et  al.  1995).  A  correlation  is  evident, 
indicating  that  the  thermal  material  in  the  filaments  which  emits  O  in  line 
emission  is  causing  free-free  absorption  of  the  radio  continuum. 


knot  on  the  near  side  of  the  Nebula  and  thus  supports  our 
hypothesis  that  the  flat  spectral  index  feature  is  due  to  free- 
free  absorption  by  thermal  electrons  in  the  line-emitting 
filaments. 

The  rms  variations  of  the  spectral  index  relative  to  the 
mean  are  2.4  <r  when  excluding  the  aforementioned  feature, 
which  indicates  that  there  is  real  structure  in  the  spectral 
index.  A  plot  of  the  intensity  of  the  [O  in]  emission  from  the 
front  of  the  Nebula,  i.e.,  that  with  negative  radial  velocity, 
against  the  spectral  index  (Fig.  6)  shows  that  there  is  a 
weak  correlation  between  the  two  in  the  sense  that  the 
is  flatter,  i.e.,  the  low-frequency  absorption  is  higher  at  the 
locations  where  there  is  more  [O  hi]  emission.  This  is  what 
would  be  expected  if  the  thermal  material  responsible  for 
the  [O  m]  emission  is  causing  free-free  absorption  at  74 
MHz. 

Assuming  an  electron  temperature  of  T  =  104  K  for  the 
optically  emitting  thermal  filaments,  our  measure  of  the 
free-free  optical  depth  constrains  the  emission  measure  to 
be  >800  cm-6  pc,  the  inequality  arising  because  not  all  of 
the  synchrotron  emission  will  originate  from  behind  the 
absorbing  feature.  Assuming  a  line-of-sight  depth  similar  to 
the  feature’s  width,  we  calculate  that  the  electron  density  in 
the  filaments,  ne  >  200  cm-3,  is  compatible  with  the  values 
expected  for  the  Crab’s  filaments  (Davidson  &  Fesen  1985). 
This  is  also  compatible  with  densities  derived  from  the 
depolarization  of  the  Crab  between  5  GHz  and  1.5  GHz 


(Bietenholz  &  Kronberg  1991)  and,  in  fact,  a  prominent 
depolarization  feature  can  be  seen  at  this  location. 


6.  SUMMARY 

We  have  presented  new,  high-resolution,  74  MHz  observ¬ 
ations  of  the  Crab  Nebula.  Combining  these  with  a  large 
body  of  observational  material  at  other  radio  frequencies 
up  to  5  GHz,  we  have  made  sensitive  radio  spectral  index 
maps  of  the  Crab  Nebula.  We  have  carefully  registered  our 
maps  at  different  frequencies,  and  we  have  taken  steps  to 
eliminate  the  effects  of  the  correlator  nonlinearity,  which 
will  tend  to  contaminate  data  of  bright  objects  taken  at  the 
VLA.  In  agreement  with  some  earlier  work,  we  find  that  the 
radio  spectral  index,  os,  of  the  Crab  is  quite  uniform  over  the 
Nebula  at  frequencies  above  300  MHz. 

We  also  find  the  following : 

1.  There  is  no  steepening  of  the  spectral  index  near  the 
edge  of  the  Nebula.  There  is  thus  no  evidence  for  any 
steeper  spectrum  shell  emission  around  the  Crab.  Such 
emission  is  expected  to  be  faint  since  the  Crab  is  still  in  the 
free-expansion  phase. 

2.  Any  gradual  steepening  toward  the  edge  of  the  Nebula 
is  small,  being  on  the  order  of  0.02  in  a.  This  is  much  smaller 
than  has  been  claimed  in  the  past.  Synchrotron  radiation 
losses  probably  produce  an  effect  somewhat  smaller  than 
this,  although  it  is  of  about  this  order  of  magnitude,  and 
hence  they  may  be  responsible. 

3.  The  spectral  index  of  the  filaments  is  no  different  from 
that  of  the  body  of  the  Nebula  for  frequencies  above  300 
MHz,  which  suggests  that  the  origin  of  the  relativistic  elec¬ 
trons  in  the  filaments  is  the  same  as  that  in  the  body  of  the 
Nebula. 

4.  We  confirm  the  existence  of  a  moving  radio  feature  in 
the  northwest  quadrant  of  the  Nebula.  It  corresponds  in 
position  to  the  dark  lane  noted  by  Oort  &  Walraven  (1956; 
see  also  Scargle  1969).  It  appears  to  be  moving  at  ~4000 
km  s-1. 

5.  At  frequencies  below  300  MHz,  we  see  evidence  of  the 
free-free  absorption  caused  by  the  thermal  material  in  the 
optical  line-emitting  filaments.  Estimates  of  the  electron 
density  in  the  line  emitting  filaments  derived  from  the 
amount  of  free-free  absorption  we  measure  are  compatible 
with  other  estimates  derived  from  the  optical  spectra  of  the 
filaments. 
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